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We examine the Chevallier-Polarski-Linder (CPL) parametrization, in the context of quintessence 
and barotropic dark energy models, to determine the subset of such models to which it can provide 
a good fit. The CPL parametrization gives the equation of state parameter w for the dark energy 
as a linear function of the scale factor a, namely w = Mo + w a (l — a). In the case of quintessence 
models, we find that over most of the wo, w a parameter space the CPL parametrization maps onto 
a fairly narrow form of behavior for the potential V (</>), while a one-dimensional subset of parameter 
space, for which w a = re(l + wo), with re constant, corresponds to a wide range of functional forms 
for V(<j>). For barotropic models, we show that the functional dependence of the pressure on the 
density, up to a multiplicative constant, depends only on Wi = w a + wo and not on wo and w a 
separately. Our results suggest that the CPL parametrization may not be optimal for testing either 
type of model. 


I. INTRODUCTION 

Observations of the accelerated expansion of the universe M indicate that approximately 70% of the energy 
density in the universe is in the form of a negative-pressure component, called dark energy, with the remaining 30% 
in the form of nonrelativistic matter (including both baryons and dark matter). The dark energy component can be 
characterized by its equation of state parameter, w, defined as the ratio of the dark energy pressure to its density: 

w=p/p , (1) 

where a cosmological constant, A, corresponds to w = — 1 and p = constant. While a model with a cosmological 
constant and cold dark matter (ACDM) is consistent with current observations, there are many models for dark 
energy that predict a dynamical equation of state. These include, for example, quintessence models, in which the 
dark energy arises from a time-dependent scalar field, cj> [§G3- (See Ref. [L| for a review), and barotropic models, 
in which the pressure is simply a prescribed function of the density fl4l - [24l |. 

In practice, an enormous number of dynamical models have been proposed for the dark energy. Hence, it has been 
considered useful to classify such models in terms of simple parametrizations for w. The simplest possibility is to 
take w to be a constant in time. However, there are a variety of problems with the assumption of constant w. With 
the exception of w = — 1 , a constant value for w does not arise naturally in the context of most physically-motivated 
models. Further, the assumption of constant w can provide wildly inaccurate results in the case that w does, in fact, 
evolve with time (25|. 

The next level of complexity is a two-parameter model for w as a function of the scale factor, a, or equivalently, 
of the redshift, z. By far the most widely-used of such parametrizations is the Chevallier-Linder-Polarski (CPL) 
parametrization, which takes re to be a linear function of the scale factor, namely [2bl f2?f : 

w = w 0 + (1 - a)w a , (2) 

where wq and w a are constants. Other two-parameter approximations have been discussed in the literature [28l. (29}: 
for a review of these and related approaches, see Ref. [3Cj. The CPL parametrization describes fairly gradual evolution 
from a value of w = wq + w a at early times to a present-day value of w = wo- It has several advantages, not least 
the fact that it is well-behaved all of the way from a = 0 to a = 1 and, for a variety of models, it can reproduce 
the predicted observable quantities (distance or Hubble parameter as a function of redshift) to extraordinarily high 
accuracy [31]. Despite the variety of proposed parametrizations, the CPL approximation has become the de facto 
standard two-parameter description for the evolution of the dark energy equation of state parameter (e.g., SHzb- 

Because of its importance, it is worthwhile to examine the applicability of the CPL parametrization in more detail. 
While it is possible to simply treat the CPL parametrization as a convenient heuristic description of the (unknown) 
dark energy component, it is also useful to understand the mapping between physical models for dark energy and 
the CPL parametrization. It has been claimed [27] that Eq. @ provides a good fit to w(a) for a wide variety of 
quintessence models. In Ref. [f32| . it was shown that thawing quintessence models with a nearly flat potential all 
converge toward the behavior given by Eq. with w a ~ —1.5(1 + wq). A similar result using different techniques 
was derived in Refs. fU, j3jj|, with w a given by w a ~ —1.58(1 + wq). This relation is consistent with results obtained 
using Monte Carlo sampling of a variety of different quintessence potentials [34| . A linear relation between w a and 
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1 + wq is also evident in the much earlier work of Kallosh et al. [H] for the special case of the linear potential. On 
the other hand, the evolution of w(a) for evolution near a local maximum of the potential (hilltop quintessence) can 
deviate strongly from a linear form |36(. 

Hence, it is interesting to determine precisely what physical models are well-described by the CPL parametrization, 
and which are not. This is, however, a difficult proposition. There are an infinite number of possible quintessence 
potentials, and even if one could sample the entire space of possibilities, the question remains of defining a “good” 
fit to Eq. ©. We therefore take the opposite approach here: we begin with Eq. © and determine, for two types of 
models (quintessence and barotropic fluids) precisely which models this parametrization maps onto. This approach 
has the advantage that one can, in a systematic way, scan over all possible values of wq and w a , and then look for 
patterns in the corresponding physical models. This mapping from the CPL parametrization to quintessence models 
has been previously explored by Padmanabhan and Choudhury [37j and by Guo, Ohta, and Zhang J3f|, both of which 
examined a single pair of wo,w a values. Barboza et al. j3£j investigated a more extensive parameter range, but for a 
different parametrization. Here we extend this earlier work by scanning over the full range of wo and w a values. We 
then perform a similar mapping for barotropic models, which have been comparatively less well explored. 

In the next section we examine the mapping of the CPL parametrization onto quintessence models and determine 
the functional form for V (</>) for this models. In Sec. Ill, we perform a similar analysis for barotropic models and find 
the pressure as a function of density for these models. Our conclusions are summarized in Sec. IV. 


II. QUINTESSENCE 


In this section, we will consider models in which the dark energy is provided by a minimally-coupled scalar field, 
</>, with equation of motion given by 


(j) + 3 He/) ——— — 0, 
d(j) 

where the Hubble parameter H is given by 

i/ = (^)=V^73, 


( 3 ) 

( 4 ) 


and the dot denotes the derivative with respect to time. Here a is the scale factor (taken to be a = 1 at the present), 
Pt is the total density, and we work in units for which 8nG = 1. Since we are interested in the evolution of dark 
energy at relatively late times, we will consider only the contributions of nonrelativistic matter (baryons plus dark 
matter) to pr, and ignore, e.g., radiation. 

The pressure and density of the scalar field are given by 

P<t‘ = Y -V(4>), ( 5 ) 

and 

P<t> = ( 6 ) 

respectively, and the equation of state parameter, w, is given by equation Cl- 

Our job is to begin with the CPL parametrization (Eq. [2]) and work backwards to derive V((f>). Such a derivation 
is provided in Refs. [33-139} (see also the earlier related discussion in Refs. [40, 4lj), so we will simply reproduce their 
results here. First, we note that Eqs. ©, ©, and © can be combined to yield 

V = \( l ~w)P4>- ( 7 ) 

Further, the density of a perfect fluid with equation of state parameter w evolves as 

-^ = -3(l + u,). (8) 

p da 

Taking w(a) to be described by the CPL parametrization (Eq. [2]), the density evolves as 

p = p < j,oi~ 3( - 1+Wo+Wa ' > e 3Wa( ' a ~ 1 \ (9) 
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where ppo is the present-day density. (The 0 subscript will refer to present-day quantities throughout). Note that 
Eqs. © and © are not specific to scalar field models; we will use them again when we discuss barotropic fluids in 
the next section. 

Combining Eqs. © and m gives an expression for the scalar field potential as a function of the scale factor: 

V(a) = ip 0O (l -w 0 -w a + w a a)a~ 3 ^ 1+Wo+Wa ^e 3Wa ^ a ~ 1 ^. (10) 

In order to express V as a function of <f>, we need an expression for (j>{a). From Eqs. ([Tj). (joj). and ©, we have 


d(f> 

dt 


^/(l + w)p<l>, 


( 11 ) 


which can be combined with the definition of H (Eq. © to give 

dcj) _ sj(l + w)p<f, 

da aH 

Substituting the CPL expression for w into Eq. (fl^l) and using H = i J(pm + /0i>)/3, with pM = Pmou~ 3 , we get 


( 12 ) 


+ Wo + Wa - W a a) 


da 


V 1 + (Pmo/ p < j,o)a 3 (' w °+ Wa )e 3wa ( 1 ~ a '> a 


(13) 


Together, Eqs. m and (fl3l) give a parametric expression for P(</>), which is our desired result. Note that Eq. (fl3l) 
yields an arbitrary additive constant in the expression for <p(a). However, this simply amounts to a shift in the value of 
(f> by an arbitrary constant in the expression for V((/>), which has no physical significance. We will take pmo/P po ~ 3/7 
in what follows. 

Before exploring the functional form for V(<p) for given values of wq and w a , we first consider the allowed ranges 
for these two quantities. Limits derived from observational data have tended to converge on a fairly narrow ellipse in 
the wq , w a plane, with negative slope (i.e., smaller wq corresponds to larger w a ) @, &]• However, the validity of these 
limits is itself dependent on the assumption that the CPL parametrization is a good fit to the evolution of w. Since 
we are interested in probing the extent to which the CPL parametrization describes specific dark energy models, we 
will allow wq and w a to vary outside of these observational limits. 

A more fundamental limit comes from Eqs. 0, ©, and ©, which together with the assumption that tp 2 /2 and 
V(4>) are both nonnegative, give the well-known bounds 


-1<w<1. (14) 

A less stringent limit comes from the requirement that the evolution correspond to a scalar field rolling downhill in 
the potential. While upward-rolling fields are physically possible, this behavior is normally a transient phenomenon 
(unless one is dealing with oscillatory behavior near the minimum of the potential [42| ). Hence, a scalar field rolling 
uphill in a potential today would correspond to yet another coincidence problem to add to the well-known approximate 
equality of matter and dark energy densities at the present. The requirement that the field be rolling downhall gives 
the constraint [43] 


a— > -3(1 — u>)(l + w). (15) 

da 

Substituting the CPL parametrization into this equation, and taking a = 1 to correspond to the present, we see that 
this translates into the limit 


W a < 3(1 - Wq), (16) 

in order that the field not be rolling uphill today. This limit is not a stringent bound, like Eq. on, so we will explore 
the parameter range outside of it, but we note that such models seem somewhat less plausible for the reasons outlined 
above. 

We now determine V{(p) for a range of values for wo and w a ■ We normalize P(</>) to its present day value, Vo, and 
we use the freedom to shift cp by an arbitrary constant to take (f> = 0 to correspond to the present. Because the CPL 
parametrization represents a heuristic fit to the data, we will not require it to hold at arbitrarily early epochs, but 
only between some initial scale factor aj and the present. The choice for is somewhat arbitrary, but we will take 









4 



-2 - 1.5 -1 - 0.5 0 

0 


FIG. 1: For quintessence models described by the CPL parametrization, the scalar field potential V((f >) normalized to its present 
day value, Vo, is plotted as a function of 4>, for wo = —0.9 (dashed, red), and w a = (top to bottom) 0, 0.5, 1, 1.5, 1.9; wo = —0.6 
(dotted, blue) and w a = (top to bottom) —0.5, 0, 0.5, 1.0, 1.5; wo = —0.3 (solid, green) and w a = (top to bottom) —1.0, —0.5, 
0, 0.5, and 1.0, where we restrict the scale factor a to lie in the range 0.3 < a < 1, and <j> is translated such that its present-day 
value is (j>o = 0. 


a,i = 0.3 in what follows, so that all of the current supernova data lie in the interval a* < a < 1. The results for V(<p) 
are shown in Figs. 1-2. Note that for the particular case of w a = 0 (so that w = wq), there is an exact solution to 
Eqs. (fTOD and (THU) , namely V ((/>) is a power of a hyperbolic sine [4J, |45|] . 

It is clear from Figs. 1-2 that the potentials corresponding to the CPL parametrization are all quite similar to each 
other. In particular, fixing wq and varying w a produces a family of functions that are nearly indistinguishable, and 
the freezing [w a > 0) and thawing {w a < 0) models correspond to potentials which appear very similar. (Note that 
some of the curves with wq = —0.9 violate the bound in Eq. (flbl) and correspond to scalar fields rolling uphill at the 
present; this can be seen in Fig. 2.) 

All of the potentials in Figs. 1-2 appear to be convex functions (d 2 V/d(f> 2 > 0). We can test this hypothesis using 
the direct expression for d 2 V/dcj) 2 from Ref. [46j] (see also Ref. 0), namely 


1 d 2 V 


(2 + 3w + q/ 2) 


1 


w 


2 


w" 3 
2(1 + w) + 4 1 


w)( 5 + 3ui + 2 q). 


H 2 d(j) 2 


1 + w 4 


1 + w 


(17) 
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FIG. 2: Expanded view of the region of Fig. 1 near <j> = 0. 


Eq. (1171) corrects a typo in the first term of Eq. (46) in Ref. [46j (E.V. Linder, private communication). Here the 
prime denotes the derivative with respect to In a, and q is the deceleration parameter, which, for a universe containing 
matter and a quintessence component with equation of state parameter given by Eq. ©, is 


13/_ w 0 + (1 - a)w a _ 

2 + 2 \ l + (/)Mo/wo)a 3( "’ 0+ ”“ ) e to “< 1 -“) 


(18) 


Taking 0.3 < a < 1, we scan over wq and w a and use Eqs. (fl7l) - (lT8l) to determine the sign of cPV/dcj) 2 . We find 
that d 2 V/d(f> 2 > 0 except for two regimes. For small values of a, d 2 V/dcj> 2 < 0 for w a > 1.5, which is outside of 
the parameter range explored in this paper. The other case for which d 2 V/d(f> 2 < 0 is for wq < —0.9, w a small and 
negative, and a near 1. We will examine this regime in more detail shortly. 

Figs. 1-2 show that the CPL parametrization maps onto a very narrow range of functions V((f>). However, this 
result seems paradoxical, since previous work [3lH33l | suggests that a wide range of functional forms for V (</>) are 
well-approximated by the CPL parametrization. The key to resolving this apparent contradiction is to note that 
these previous studies did not indicate that quintessence models map to the full wq, w a plane. Instead, Refs. fdTI [dd| 
show that thawing quintessence maps to a one-dimensional subspace of this plane, namely, the line defined by the 
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relation 

w a = k(1 + w 0 ), (19) 

with n variously estimated to be in the range from —1.6 to —1.5. 

To understand this result, we assume the validity of Eq. (1191) and derive the corresponding forms for V(4>) that 
satisfy it. These forms for V(<f>) are shown in Fig. 3 for k ranging from —0.5 to —3. Note that in this case, we 
choose a, to be the value of a for which w(a) = —1; this is the earliest allowed initial value for a for which this 
parametrization can be valid. This yields a* = 1 + 1 /k, independent of wq and w a - 



FIG. 3: As Fig. I, for w a = k(1 + wo), with wo = —0.9 and k =(bottom to top) —0.5 (green), —l (blue), —1.5 (red), —2.0 
(black), —2.5 (cyan), —3.0 (magenta), where <j> is translated such that its present-day value is (j >o = 0, and at is the smallest 
value of a for which w(a ) > — 1. 

These forms for V((f>) become nearly linear for n = —1.5, while diverging from linear behavior on either side of this 
value. This suggests that the particular choice of Eq. (usd corresponds to a nearly linear potential. Further, when w 
is close to —1, (which is always the case when Wo « — 1 and w a < 0), (f> rolls only a short distance down the potential, 
and any sufficiently smooth potential will appear to be linear on small enough scales. This explains why thawing 
potentials with w near —1 are described so well by Eq. (USD- 

Further, Fig. 3 corresponds precisely to the potentials for which V(<j>) is concave {(PV/dcf 2 < 0) for a near 1. This 
corresponds to <f> near 0 in Fig. 3. It is clear from this figure that d 2 V/d(j > 2 is negligibly small in these cases. Thus, 
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we can conclude that for the parameter ranges examined in this paper, d 2 V/d(j) 2 is never large and negative. Thus, it 
makes sense that w(a) for the hilltop models examined in Ref. [36| (for which d 2 V/d<p 2 is large and negative) strongly 
diverges from linear evolution in a. This does not mean that the CPL parametrization is a poor fit for all models with 
d 2 V/d(j) 2 < 0. If d 2 V/dq L > 2 is sufficiently small, then our argument from the previous paragraph can apply, yielding 
evolution well-described by the CPL parametrization and Eq. (fl9l) [H, (H, [36]. 

Thus, we conclude that over most of the range in wo and w a , the CPL parametrization maps onto a very narrow 
form of behavior for V (</>), while a narrow range in the wq, w a parameter space (given by Eq. 1191) maps to a wide 
range of quintessence models, namely thawing quintessence with w near —1. However, it is precisely these models 
that appear to be allowed by current observations, as emphasized by Ref. [33]. 


III. BAROTROPIC MODELS 

In this section, we examine barotropic models, for which the pressure is a fixed function of the density: 

P = f(p)- (20) 

Particular models of this form include the Chaplygin gas ai and the generalized Chaplygin gas [b6j, the linear 
equation of state [l7], EH and the affine equation of state Hi (note these are actually the same model), the 
quadratic equation of state [l9| , and the Van der Waals equation of state [2l|, [22| . A general study of the properties 
of barotropic models for dark energy was undertaken in Ref. [23} and further extended in Ref. [53 . Note that there 
is a simple mapping between the barotropic models discussed here and kinetic k-e ssence models [231 ] . so the results 
presented here can be extended in a straightforward way to the latter set of models. 

Ref. [23| provided two reasonable constraints on the behavior of w(a) in barotropic models. To prevent instabilities, 
the sound speed cl = dp/dp must obey cl > 0, while causality requires cl < 1. While each of these limits translates 
into a constraint on w(a), it is easier in this case to derive the form for p as a function of p and then apply the direct 
limits on dp/dp: 


0 


< 


dp 

dp 


< 1 . 


( 21 ) 


If the equation of state for the barotropic fluid is given by Eq. ©, then the expression for the density p given by 
Eq. © is valid as well, i.e., 


p = poa -3(i+^o+v, a ) e 3 Wa (a-i) ; 


( 22 ) 


Further, the pressure will be given by wp, so 

p = Po[w 0 + (1 - a)w a }a- 3 ^ +W0+w ^e 3w ^ a - 1) 


(23) 


Eqs. (l22l) and (l23l) together provide an implicit expression for f(p) in Eq. (l20l) . As in the previous section, we will 
not require the CPL parametrization to hold all of the back to a = 0, but only for etj < a < 1, where we take etj = 0.3. 
In Figs. 4-5, we plot p/po as a function of p/po for the scale factor lying in this interval. 

The behavior of p as a function of p for these barotropic models obeying the CPL parametrization displays a much 
wider set of functional behaviors than is the case for V((f>) in the quintessence models. Note, however, that all of the 
curves converge on the same value of p/ pn (= wq) at p/po = 1, corresponding to a = 1. 

We can apply the limits from Ref. 1231 given by Eq. m directly to the models displayed in Fig. 4. It is clear that 
causality (dp/dp < 1) is satisfied for all of the models under consideration. However, the requirement that dp/dp > 0 
is violated by some of these models. From Figs. 4 — 5, it is clear that causality requires relatively large positive values 
of w a (i.e., freezing models rather than thawing models). 

We can make a stronger statement about the allowed forms for f(p) and their dependence on wq and w a ■ Eq. © 
can be used to express the scale factor as a function of w, wq, and w a : 


w a + w 0 - w 

Wa 


(24) 


and this expression can be substituted into equation (E^l) to give 


P = 


p 0 \w a \ 3il+Wi) e~ 3wa \wi - 1i ,|-3(l+t»«) e 3(«ii-lB) j 


(25) 






10 20 30 40 50 

P / Pq 


FIG. 4: For barotropic models described by the CPL parametrization, the pressure p, given as p/po, is plotted as a function of 
the density p, given as p/po, where pa is the present-day dark energy density, for wo = —0.9 (dashed, red), and w a = (bottom 
to top) 0, 0.5, 1, 1.5, 1.9; wo = —0.6 (dotted, blue) and w a = (bottom to top) —0.5, 0, 0.5, 1.0, 1.5; wo = —0.3 (solid, green) 
and w a = (bottom to top) —1.0, —0.5, 0, 0.5, and 1.0, where we restrict the scale factor a to lie in the range 0.3 < a < 1. 


where we define Wi to be the value of w at a = 0 in the CPL approximation, namely, wi = wq + w a ■ Note that Eq. 
m is an implicit equation for p(p), since w =pjp. Now define 


P = P/PG, 

p = p/pc, 


with the constant density pc given by 

Pc = Po\w a \ 3 { 1 +Wi) e~ 3w ' 


Then Eq. (E5l) becomes 


p= | Wi-p/pl-^+^e 3 ^-^. 


(26) 

(27) 


(28) 


(29) 


From this equation, it is clear that the expression for p as a function of p depends only on Wi and not on w a and wq 
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independently. Writing this function as p = g(g 5), we then have 

p/pc = g(p/pc) (30) 

where the function g depends only on Wi. Thus, aside from an overall multiplicative constant in both p and p (which 
is the same for p and p, so that w is unaltered), f(p) in Eq. (EU1) depends only on Wj, and not on w a and wq separately. 


IV. CONCLUSIONS 

Lacking a priori knowledge of the actual physical model underlying the accelerated expansion of the universe, it 
is impossible to make blanket statements about the utility of of the CPL parametrization for approximating such 
models. However, we can make statements about individual dark energy models; the implicit question is the extent 
to which a linear parametrization of w provides useful information about the parameters of a given physical model. 

For the case of quintessence, we have seen that the full parameter space of the CPL parametrization maps to only 
a fairly narrow form of functional behavior for V((/>). However, a one-dimensional subset of CPL parameter space, 
namely w a = k(1 + wq), with k ps —1.6 —1.5 corresponds to a wide range of functional forms for V(cf>). This issue has 
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been noted by Linder [33j , who pointed out that such models are better fit by the one-parameter model of Eq. m 
than by the full CPL parametrization. Further, Ref. introduced a more accurate two-parameter approximation, 
using Eq. (flT)l) as a starting point, that better fits the behavior of thawing quintessence models. Our results provide 
further evidence in favor of this approach. 

We note a similar situation for barotropic models. For these models, the CPL parametrization does a much better 
job of capturing the full range of possible functional forms for p = f(p ) (the pressure as a function of density). 
However, we have shown explicitly that the barotropic models that correspond exactly to the CPL parametrization 
belong to a one-dimensional subset of the full CPL parameter space; namely, they are fully determined (up to a 
multiplicative constant) by the value of Wi = wq + w a - Thus, the full CPL approximation may not be the most useful 
way to characterize such models. 

Considered as a straightforward linear approximation for the equation of state parameter as a function of the 
expansion factor, the CPL parametrization is undeniably useful, particularly when attempting to detect deviations 
from the ACDM model. However, as a probe of physical models which might correspond to non-ACDM behavior (or 
at least the types of physical models considered here), the CPL parametrization is not optimal, and, at least for the 
case of quintessence, it seems more appropriate to move toward the kind of parametrization discussed in Ref. [33]. It 
would, of course, be useful to extend this study to some of the other classes of models for dark energy that have been 
proposed in the literature. 
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